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Q ' Abstract 

The calculation of the third order susceptibility still is a long standing fundamen- 

> 

' tal problem of particular importance in nonlinear nanooptics: Indeed, cancellation 

00 : 

. of size-dependent terms coming from uncorrelated excitations is expected, but up to 

now shown for very simple Hamiltonians only. Using a many-body theory recently 

m : 

. developed to handle interacting close-to-bosons, we prove it here for arbitrary H. 

I This new formalism actually provides the first clean way to calculate nonlinear 

susceptibilities, with results different from previous ones. 

o 
o 

> 
X 

d ' PACS.: 71.35.-y Excitons and related phenomena 
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Size enhancement effect in nonlinear optical processes is of great interest in nanoscience 
and technology as it provides an additional degree of freedom for the control of matter. 
Care is however necessary in the calculation of nonlinear susceptibilities due to possible 
cancellation among their various terms [1] . Indeed, for non-interacting bosons, the cancel- 
lation is so complete that all nonlinear susceptibihties vanish identically. This shows that 
optical nonlinearities only come from the non-bosonic nature of the excitations, and/or 
their possible interactions. 

In the various terms of the third order susceptibility x^^^ ■, the ones with uncorrelated 
excitations increase linearly with the sample "volume" . As the susceptibihty is an 
intensive quantity, all these volume-linear terms have clearly to cancel. However, along 
with them, other terms in L° can also disappear. This is why, if this cancellation is not 
handled properly, the obtained x*-^-* may well be incorrect. The cancellation problem is 
thus crucial for nonlinear optics [2]. 

A lot of controversial arguments [3] have been raised to justify the cancellation of 
volume-linear terms. Some fundamental aspects were clarified through the study of a 
simple model consisting of non-interacting ID Frenkel excitons with Pauli exclusion. It 
allows to analytically show [4] the cancellation of these volume-linear terms for arbitrary 
frequencies. Numerical evaluation of the remaining terms shows a size enhancement for 
sample sizes within the long wavelength approximation and a possible saturation when 
the size approaches the smaller of either, the coherence length due to non-radiative scat- 
terings, or the relevant wavelength of light [5]. This model calculation, extended to higher 
dimension for excitons having some kind of interactions [6], again shows cancellation of 
the volume-linear terms. 

Although these simple models give a reasonable picture of the problem, they are not 
sufficient to carry out reliable calculations in realistic situations. A general method for 
calculating nonlinear susceptibilities which does not suffer from this cancellation problem, 
was thus highly needed. 

A new many-body theory has been recently developed [7-13] to handle interactions 
between close-to-bosons: It allows to extract from the quantities of physical interest, the 
parts coming from pure Pauli "interaction" and the parts coming from Coulomb interac- 
tion, either through direct or exchange processes. We use this "commutation technique" 
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here, to solve the cancellation problem of the third order susceptibility, by proving that 
the terms independent from exciton interactions cancel out exactly for any matter Hamil- 
tonian. As expected, the remaining terms of x^^^ depend on interactions only, pure Pauli 
or Coulomb dressed by exchange, our new formalism allowing to write them in a compact 
way. 

The third order nonlinear susceptibility is the integral kernel relating the field vector 
potentials to the current density induced by the semiconductor-photon coupling taken at 
third order. For photons turned on adiabatically from — oo to to, this induced current 
density reads [14] 

J = r dt, r dt2 r dts{v\[[[io,hii2],i3]\v) , (i) 

J —oo J —oo J —oo 

where Ij^o — —ie^^*^Hl^^e~''^^^, while Iq = e*^*° Jre~*^*°. H is the matter Hamiltonian, 
Jr the induced current density operator and Hl"^ — e^*(e*'^*C/ + h.c.) the couphng between 
matter and (a;, Q) photons, if we consider one photon field only for simplicity, e = 0+ 
being the adiabatic switching factor. Jr and H"^^ create or destroy one excitation in the 
sample. In terms of creation operators for the one-excitation eigenstates Bj^ defined as 
(H — Ejn)Bl^\v) — 0, they formally read 

Jr — Jr,m ) Jr,m — jm(j^) -^m + h.C. , (2) 

m 

U C^Yl l^mBm , (3) 
m 

eq. (3) being valid within the rotating wave approximation. For Wannier excitons, m 
stands for {I'mjQm) with being the exciton center of mass momentum and z/^ the 
relative motion index; //^ = A.G5q^_q (r — Olum) L^^'^, while imi^) — — Ge'*^*" "" (r = 
Olvrn) L'^/^, with A being the field vector potential and G the valence-conduction Kane 
vector [15]. 

Due to eq. (2), J thus appears as Z^mJm, where is given by eq. (1), with Iq 
calculated using Jr,™., instead of Jr. By expanding the commutator of eq. (1) and by 
noting that jj_^Q = —Ij while = lo, Jm actually reads Jm + c.c, with J^ = K0123 + 
K2103 + K3102 + K3201 and 

Kjjfei ^ [° dti [' dt2 dts {v\Ii Ij Ik Ii\v) . (4) 
J —00 J —00 J —00 
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The purpose of this letter is to show how the 8 terms of Jm can be calculated formally, 
i. e., without knowing H, in order (i) to prove that the sum of their non-interacting 
contributions cancel out exactly, (ii) to write a compact expression of the remaining terms. 

As exphcitly shown below, the /j //|v)'s of eq. (4) split into "easy" terms in 
which the middle state is the vacuum — so that they can be calculated exactly as they 
depend on one-excitations only — and "tricky" terms in which the middle state contains 
two excitations. As the two-excitation eigenstates are usually unknown, it has long seemed 
hopeless to prove that, whatever H is, these "tricky" terms contain parts which exactly 
cancel the "easy" terms. 

Using the "commutation technique" [8,9] we recently developed to handle many-body 
effects between close-to-bosons, we can perform commutations between Bl^ and the in- 
dependent operators of these "tricky" terms, to rewrite them as a sum of "non-interacting" , 
"Pauh" and "Coulomb" contributions. The "tricky" terms of the Kyfci^oi23's of then 
read j;,(r, to) (Pijki + lijki + Sijki), with j;,(r, to) = j^(r)e(^'+''^)*o, the prefactors f^ijM, lijkh 
Sijki coming from non-interacting, Pauli and Coulomb processes, respectively. Kq/^s'^-^'* 
has a similar form, with jj!„(r, to) replaced by jm (1*5^0)) while the "easy" terms of all the 
Kijki's of Jm read j„(r, to) aijki- 

Prom the formal expressions of these a's and /?'s, which appear as H matrix elements, 
it is then possible to show the far from obvious relations, ci;oi23 + Q!2io3 + (^tio2 — 0) 
^3201 + ^3102 + Plmi = 3.nd /9oi23 + /?32oi = 0) whatever H is [16]. Consequently, only 
remains in the induced current density, the 7 and 5 parts coming from Pauli and Coulomb 
interactions between two excitations, as expected. It ultimately reads 



- JmliJ e ^1 ^, ^, 



X 



m p,q,n q ■ 

edir tin 

\ . ipqmn ipqmn , 



+ C.C. , (5) 



El = Ei — CO being the detuned exciton energy, the next term having one more ^ divided 
by detuning. 

\qmn, Cpgmn C^mn exchange, direct Coulomb and "in" exchange Coulomb 

scatterings of the "commutation technique" . Defined through [8,9] 

[Bm, Bj] = Smi - Dmi, [Dmi, -^j] = 2 ^ ^mnijBl, (6) 
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— - . (11) 

-cx) 2m X — H — in 



[H, Bj] = eM + y/, [V^;, S]] = ^ ^^ElEl, (7) 

mn 

and i^^nij — Y.rs ^mnrs^tfij-i they read in terms of the one-excitation eigenwavefunctions of 
H (see eqs. (26,28) of ref. [8]). The "deviation-from-boson" operator Dmi and the Xmnii^ 
physically come from Pauli exclusion which makes the excitons close-to-bosons only, while 
the "creation potential" and the ^mmj's come from Coulomb interactions with the i 
exciton. 

To calculate J^, we have also used [10] 

which follows from eq. (7) and 

e-'"^Bl = Bl e-(^+^"')- + wI{t) , (9) 

/+00 p-'(:r-iri)T i 
— — = , (10) 
-oo 2i7r x-H-iri"' x-H- Em -iv 

which follows from eq. (8) and the integral representation of e"*^"^ which, for r < and 
?7 > 0, reads [17] 

IT] 

To show how this "commutation technique" works and to grasp the physics which 
controls the a, /?, 7, 5 prefactors, let us consider a 'i^ijki particularly complicated to 
calculate, namely K2103, because in addition to the fact that its middle state has two 
excitations, the times tj are "shaken up" compared to K0123, which makes the integration 
over times more difficult. 

e^*^*|v) = \v) and U\v) = 0, so that {v\l2 and /aju) reduce to one term. As the non- 
zero contributions of IiIq to k2io3 = {v\l2lilQh\v) are the ones in U^Bm and C/-B^ only, 
k2io3 splits in two. Its "easy" term, containing U^Bm, reads k2ioT'' = i e^^^'^i+^'^+'^^ 02103 
jm(r,to) with 

02103 = {v\Ue-'"'^' e-^^^i B^ e^H' (n+r^+ra) ^t|^) ^ (12) 

where H' — H — cu, while the Tj's, defined ets ti — Ti + ti-i, run from —00 to 0. The 02103 
middle state being the vacuum, the middle H can be replaced by and the right H' by 
E'^. As {v\BmW\v) = fi*^, the formal integration of this "easy" term over the r^'s readily 
gives a Q!2io3 jm(r, ^o) contribution to K2103, with 

q;2103 = 7 ^ {v\U i U^v) . (13) 
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The k2io3 "tricky" term reads k^\of^^ = i e<3ri+2T2+T3) 62103 j;,(r, to) with 

{v\U e-'""'^'' U e-'^""^ Bl e^^'(^i+^2+r3) ^t|^^ ^ (14) 



^2103 



and H" = H — 2u. The integrations over ri and T2 are not straightforward because these 
times appear in different places. We can put the Ti's together by passing e~*^"^i over Bj^ 
through eq. (9). This sphts 62103 as C2103 + <i2io3- From the first term of eq. (9) we get 

C2103 = e-'""'^^' {v\Ue-''''^' U Bl e''''^^'+^'^ U^v) . (15) 

T2 still is in two places. To go further, we rewrite UB^ as BlU + [U, Bj^], i. e., BIJJ + 
A*m — J2qfJ'qDqm, duc to cqs. (3,6). In the part with BlU, we can replace the left H' by 
E'^ since the middle state is then the vacuum. As {v\UBHv) — the two first terms 
of UBl give a /?2io3 *o) contribution to K2103, with 

P2IO3 — 



E' + 3ie 



x{v\u( ^- -] \v) . (16) 

As for the third term of UB^, in Dg^, it physically comes from Pauli "interaction", so 
that it is zero for boson-excitons. The contribution it induces to K2103 reads 72103 jml'^) ^0) 
with 



00 

iH'T2 



x{v\U e-'"'-' D,m -^rz— 1^) • (1^) 



Since H' acts on one-excitations only, the integration over T2 can be performed exactly, 
using eq. (3). As {v\BpDqmBl^\v) — 2\pqmn, due to eq. (6), we find [18] 

72103 - ^ z.^ _ 2^^) • I ; 

In passing e~*^"'^i over Bj^ in eq. (14), we also generate a 0^2103 term which comes from 
W^{ti) in eq. (9), i. e., Coulomb interactions with the m exciton. This term produces a 
52io3jm('')^o) contribution to K2103 with 



-'2103 



i I dT2 e^'^' {v\ U e-'^'^^ U Wl \v) , (19) 

J —00 11 — Id 



where W^^ is linked to VF^(ri) through 

Wl = -i f dn Wl(ri) e(3^+^(^'+2'-))^i . (20) 

J—oo 

Using eq. (10), it reads 

J-oc 2i7ry-H"-ir] 

X - (21) 

{y-H' + ie'){y-H'-E'^-irj) ' ^ ^ 

with ?7 chosen such that e' = 3e — 77 > 0, to insure convergence for ti —00. 

The integration over T2 can be performed exactly, using again eq. (3). This leads to 

[18] 

r _ Wpqran l^p (^q t^n ('221 

~ E'^ + 3te {E'^ - ze) {E^ - E'^ - 2ie) " ^ ^ 

Wpqmn — {E'^ + 3ze) {v\BpBqWl^Bl^\v) corresponds to Coulomb processes between the m 
exciton and the photocreated ones {p, q) and n. When inserted in Wpqmn, the wl^, of eq. 
(21) can have its two right H' replaced by E'^. Integration over y, done by residues, yields 

Wpqmn = (^^l^^ Bq _ ^„ _ BUv) . (23) 

Wpqmn Cannot be calculated exactly because H" acts on two excitons. It can however be 
expanded in Coulomb interactions using eq. (8). At lowest order, this leads to replace H" 
by Ep + Elj. Prom eq. (7) and the scalar product of two-exciton states, {v\BpBqBlBl\v) — 
5pr5qs + 5ps^qr — ^Xpqrs, which foUows from cq. (6), we end with 



tdir tin 

Wpqmn ' ^ E'^ - E'^ - E'^ - Zte ^ ' 



Cpqmn exchange Coulomb scattering with its Coulomb interactions between the 
"in" excitons {m,n) only. Once again [8,9,10], it differs from the "out" exchange Coulomb 

scattering ^°q^n — Yjr,s ^pqrs -^rsmn appearing in the effective bosonic Hamiltonian for 
excitons [19], widely used up to now, in spite of the fact that it is unphysical because not 
hermitian. 

The "commutation technique" being new, and the project to extract from J its non- 
interacting terms ambitious, it seemed to us useful to detail the calculation of a particular 

Jiijki, even if this calculation may appear somewhat technical. The expressions of a, (3, 
7, 5's it provides, are actually quite interesting. 



(i) The a's and /?'s, independent from Pauli and Coulomb interactions, contain three 
exciton-photon coupUngs fj,pfj,qii^. In the case of Wannier excitons, they generate three 
factors L^/^ and force Qp = = Qn = Q so that no sample volume comes from the 
sums over (p, n). By including the L"^/'^ factor of j^(r), these cu's and /3's thus generate 
volume- linear contributions to J^, in agreement with model Hamiltonians. They have 
to and do cancel exactly for any kind of excitons. 

(ii) The 7's and 5's have the same number of exciton-photon couplings. They however 
have one additional Pauh "scattering" for 7's, and Coulomb scattering for 5's. As, for 
Wannier excitons [8,9,12], these scatterings both behave as the exciton volume divided by 
L^, the 7 and 5 contributions to the current density are indeed sample volume free. 

(iii) Finally, at large detuning, the induced current density is dominated by pure Pauli 
processes: As obvious from eqs. (6,7), the Xpqmn^ are dimensionless, while the ^pgmn's 
have the dimension of an energy, so that the (5's do have one more energy denominator 
than the 7's, which makes them smaller in the large detuning limit. 

Conclusion 

Using a new many-body theory for interacting close-to-bosons, we have succeeded in 
extracting the various parts of the induced current density linked to uncorrelated excita- 
tions and proved their cancellation, without knowing H. This formalism actually provides 
the first clean way to calculate nonlinear susceptibilities. After combining all its con- 
tributions, the induced current density appears in terms of the exchange and Coulomb 
"scatterings" of the "commutation technique" (see eq. (5)). Explicit calculations for Wan- 
nier and Frenkel excitons, as well as polarization effects linked to spin degrees of freedom, 
will be presented elsewhere. We can however say, just from dimensional arguments, that, 
at large detuning, nonlinear susceptiblities are entirely controlled by Pauli interactions 
between excitations, without the help of any Coulomb process — result beyond the reach 
of the effective bosonic Hamiltonians for excitons. 

This work was supported by the CNRS of Prance and by the Grant-in- Aid for Scientific 
Research (No. 15540311) of the Ministery of Education, Culture, Sports, Science and 
Technology of Japan. 
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